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Abstract 



The development of accurate theoretical failure, fatigue, and wear models for ultra-high molecular 
weight polyethylene (UHMWPE) is an important step towards better understanding the 
micromechanisms of the surface damage that occur in load bearing orthopaedic components and 
improving the lifetime of joint arthoplasties. Previous attempts to analytically predict the clinically 
observed damage, wear and fatigue failure modes have met limited success due to the complicated 
interaction between microstructure deformations and continuum level stresses. In this work, we 
examined monotonic uniaxial and multiaxial loading to failure of UHMWPE using eight failure 
criteria (maximum principal stress, Mises stress, Tresca stress, hydrostatic stress, Coulomb stress, 
maximum principal strain, Mises strain, and chain stretch). The quality of the predictions of the 
different models was assessed by comparing uniaxial tension and small punch test data at different 
rates with the failure model predictions. The experimental data was obtained for two conventional 
(unirradiated and gamma radiation sterilized in nitrogen) and two highly crosslinked (150 kGy, 
remelted and annealed) UHMWPE materials. Of the different failures models examined, the chain 
stretch failure model was found to most accurately capture uniaxial and multiaxial failure data for both 
the conventional and the highly crosslinked UHMWPE materials. In addition, the chain stretch failure 
criterion can readily be calculated for contemporary UHMWPE materials based on available uniaxial 
tension data. These results lay the foundation for future developments of damage and wear models 
capable of predicting multiaxial failure under cyclic loading conditions. 



Key Words — constitutive modeling, failure modeling, ultra-high molecular weight polyethylene, 
damage, multiaxial mechanical behavior. 
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1. Introduction 



Mechanical failure of ultra-high molecular weight polyethylene (UHMWPE) components is 
one of the major factors controlling the lifetime of total joint arthoplasties. The failures that occur in 
vivo are often introduced by contact damage in the form of wear and fatigue [18] caused by relative 
sliding motion between different components, primarily due to metal-to-polymer interactions. To 
increase the long-term performance of total joint arthroplasties, investigators have attempted to 
understand the micromechanisms underlying surface damage in UHMWPE components [29, 30]. The 
reasoning for this approach is that a fundamental understanding of the micromechanisms of surface 
damage will facilitate predictions of performance on the component level, which in turn can enable 
design and material improvements. Despite these efforts, the micromechanics controlling the 
evolution of surface damage in UHMWPE components are still not fully understood and the state of 
the art in modeling the failure events is far from complete. 

Modeling the damage and failure of UHMWPE, based on micromechanical events, is a 
challenging task. On the local scale, the material consists of intermixed amorphous and 
semicrystalline domains. The material behavior is further complicated by the evolving 
micromechanical state due to applied external loading. Failure of UHMWPE on the microscale is 
driven by molecular chain scission [12, 17]. The progressive failure events are promoted by the 
magnified stress and strain state in the neighborhood of already broken molecular chains, resulting in 
damage localization that is manifested by growth and coalescence of microvoids (also referred to as 
cavitation) [25, 31]. 

In general, failure of any polymeric component can be defined as the inability to perform an 
intended function. Examples of potential failure modes are: excessive plastic deformation, static 
fracture, fatigue fracture, and wear. This paper focuses on static fracture, which is here defined as loss 
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of material continuity through the formation of new surfaces within the UHMWPE component. This 
failure mode is the most fundamental, and any failure model of common clinical failure modes, such as 
fatigue or wear, will need to also address failure under monotonic loading conditions 

There are two general classes of failure models for polymers: isostate models and evolving 
state models. An isostate failure model is a failure model in which the material state is assumed not to 
change during the applied deformation. The most general isostate failure model is given by 

/(T,E) >f C rih (1) 

where T is the current stress state (expressed, for example, in terms of the Cauchy stress tensor), E is 
the current strain state (expressed, for example, in terms of the logarithmic strain tensor), /(•) is a pre- 
specified failure model, and f cr u is a material parameter specifying the critical value of the failure 
condition. There are two main classes of evolving state models: state variable and history-based 
models. A general failure condition in state variable form can be written: 

/(£,T,E)>/ cm , (2) 

where £ ( i e [1,/V]) is a set of state variables specifying the current state. The state variables are 
assumed to evolve with the deformation according to the rate-equation: 4= 8 (^ ,T,E) , where g(-) is a 
pre-specified evolution equation. 

An evolving state model in history-based form can be written 

\f(T(t’),E(t’))dt’>f mt , (3) 

where Equation (3) symbolizes that the failure not only depends on the current values of stress and 
strain, but also on their complete time histories. One interpretation of the history-based integral model 
can be obtained by defining the integral on the left-hand side of Equation (3) as the damage D(t). The 
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rate of change of damage is then given by D(t)= f { T,E), and the history integral representation can 
consequently be seen as a special case of the state variable-based evolving state failure model. 

The goal of this study was to examine different approaches for predicting static fracture, and 
then in a later study use this necessary information to address fatigue and wear failures. In this work, 
we have used monotonic uniaxial and small punch (equibiaxial) tests to failure to examine the relative 
capabilities of different failure models to predict fracture for two conventional (unirradiated and 
gamma radiation sterilized in nitrogen) and two highly crosslinked (150 kGy, remelted and annealed) 
UHMWPE materials. Currently, conventional UHMWPE is the dominant material for hip and knee 
total joint replacement implants. However, highly crosslinked UHMWPE materials are gaining in use, 
particularly for hip replacement components, due to improved wear resistance [19]. 



2. Materials and Methods 

There are a number of different types of models that can potentially be used to predict failure of 
UHMWPE, including: stress-based models; strain-based models; energy-based models; linear elastic 
fracture mechanics (LEFM) models; and damage models [6]. The first four types of models are 
constant state models, whereas the damage models are evolving state models. It is unlikely that an 
energy-based model will be predictive since large strain energy can be stored in hydrostatic 
compression but the material will not fracture in this mode. LEFM models are a specific case of a 
stress-based constant state model. Thus, in this study, we examined the following eight isostate failure 
models: (1) maximum principal stress; (2) Mises stress; (3) Tresca stress; (4) hydrostatic stress; (5) 
Coulomb stress; (6) maximum principal strain; (7) Mises strain; and (8) chain stretch. Most of these 
models are well known and are discussed in detail in most introductory solid mechanics textbooks 
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[20]. For completeness, three of the lesser known models (the Coulomb model, the Mises strain 
model, and the chain stretch model) are discussed herein in more detail. 

The Coulomb stress is defined as 

° Coulomb = fa — <T 3 ) + H ((7i + C7 3 ) (4) 

where o\ is the largest principal stress, <73 is the smallest principal stress, and ^= 0.10 is taken as the 
friction coefficient. The Mises strain is defined analogously to the traditional Mises stress, i.e. 

^Mises = ^[(^, - ^2 ) 2 + (£ 2 ~ £3 f + fe ~ £< f ] ■ (5) 

where e, are the principal strains. 

The chain stretch failure condition is motivated by our recently developed physically-inspired 
constitutive model for UHWMPE, the Hybrid model (HM) [9, 10]. The HM is based, in part, on some 
of the constitutive modeling work by Boyce et al. [2, 3, 16] who have successfully modeled, glassy 
amorphous polymers, semicrystalline polymers, and rubbers. However, these theories did not 
adequately describe UHMWPE; thus, the HM was developed and validated specifically for 
semicrystalline UHMWPE. One of the key features of the HM is the concept of molecular chain 
stretch [7]. In the chain stretch failure model, an estimation of the average molecular chain stretch is 
given by the expression 

Khain = ^\[ e2£l +e2£l +<?2£3 ] ’ (6) 

where, as before, r.\, £ 2 , and £3 are the principal strains. This representation enables a mapping from the 
applied deformation state to the local molecular stretch state, and has been examined using molecular 
dynamics and Monte Carlo simulations [ 8 ]. In the chain stretch model, failure is assumed to occur 
when the average molecular chain stretch reaches a critical level. 
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To evaluate the quality of the predictions of the different failure models, it is necessary to use 
appropriate experimental test methods and well-characterized materials. In this study, we examined 
four different types of UHMWPE obtained from ram-extruded rods of GUR 1050 that were subjected 
to different radiation and heat treatments. The four materials [20] were: (1) as-received unirradiated; 
(2) gamma radiation sterilized in Nitrogen with 30 kGy; (3) gamma radiation treated with 100 kGy, 
giving a highly crosslinked microstructure, and subsequently annealed at 110°C for 2 hours; and (4) 
gamma radiation treated with 100 kGy, giving a highly crosslinked microstructure, and subsequently 
remelted at 150°C. The microstructures of the materials studied in this work have been extensively 
examined elsewhere [20]; e.g., the degree of crystallinity of the three materials was determined to be 
0.50 for the unirradiated, 0.51 for the sterilized material (30 kGy y-Ni), 0.61 for the crosslinked 
material heat treated at 110°C, and 0.46 for the crosslinked material heat treated at 150°C. The four 
materials were tested in uniaxial large strain monotonic tension to failure and in small punch testing to 
failure. 

2. 1 Uniaxial Tension to Failure Experiments 

The stress-strain behavior from the uniaxial tension to failure tests has been reported previously 
[16] for the four materials examined in this study. Uniaxial tension behavior of other types of 
conventional and crosslinked UHMWPE have been reported in other studies [1, 13, 15, 21-23, 27]. 
Briefly, the uniaxial tension to failure experiments used specimens that were prepared and tested 
according to ASTM D638 [4], The testing was performed using a 4204 Instron electromechanical load 
frame, equipped with a 5 kN load cell. The specimens were tested under displacement control at 
crosshead displacement rates of 30, 75, and 150 mm/mi n, approximately corresponding to initial true 
strain rates of 0.01/s, 0.02/s, and 0.05/s. A non-contacting video extensometer was used to measure 
axial strain [20]. Five specimens were tested for each experimental condition, corresponding to a total 

of 60 specimens (5 specimens per material x 4 materials x 3 rates). 
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2.2 Small Punch Testing to Failure 



The load-displacement behavior from the small punch failure tests has been reported previously 
[18]. Briefly, the small punch testing was performed on disc shaped specimens measuring 6.4 mm in 
diameter and 0.5 mm in thickness. The specimens were tested using an MTS Model 859 testing 
machine by indentation with a custom-built, hemispherical head punch at three different punch 
displacement rates (0.5 mm/min, 1.0 mm/min, and 5.0 mm/min). The tests were performed at room 
temperature (20°C ± 2°C), according to ASTM F2 183-02 [5]. Five specimens were tested for each 
experimental condition, corresponding to a total of 60 specimens (5 specimens per material x 4 
materials x 3 rates). 

The experimental data obtained from the small punch tests is in terms of force and 
displacement of the punch head. This data can be converted to stress and strain by simulating the 
small punch tests using the finite element (FE) method. To perform these FE simulations an 
axisymmetric representation of the experimental setup was created (Figure 1). In this FE 
representation, the test specimen was modeled with 360 quadratic axisymmetric elements, and the 
punch and fixture were assumed to be rigid. The FE simulations were performed using the commercial 
finite element package ABAQUS (HKS Inc., RI). 

To accurately determine the stress and strain distributions in the test specimens, it is necessary 
to have an accurate constitutive model capable of predicting the material behavior in large-strain 
multiaxial deformation states. In this work we used the HM [9-11]. The HM has been calibrated and 
validated for the four UHMWPE materials that are used in this study. The material parameters that are 
used in the simulations of the four materials and a brief summary of the HM are presented in Appendix 
A. 
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To enable failure predictions of the small punch testing, a total of 12 FE simulations were 
performed (4 materials x 3 rates). In these simulations the punch was displaced to the maximum 
displacement value that was experimentally observed in the five repeated experiments on that specific 
material and applied displacement rate. The maximum stress and strain values at the experimentally 
observed failure states for each of the five repeated tests were then extracted from the FE simulations 
by finding the location having the maximum stress and strain and extracting the stress and strain 
tensors at that location. 

2.3 Failure Modeling 

The failure parameters for each material (e.g. the Mises stress at failure) were obtained from the 
average of the experimentally observed values from all uniaxial and small punch tests. The quality of 
the different failure models was evaluated by comparing the predicted failure values for the two 
loading modes at each of the three displacement rates. The quality of the predictions was quantified 
with the coefficient of determination [9], r 2 , between the experimental failure values and the predicted 
failure condition for each material. The predicted failure condition was assumed to be independent of 
loading rate. 

The predictive capability of each model was also graphically represented using “box and 
whiskers” plots. In these plots, the center 50% percentiles of the data are represented using a box. The 
median is represented with a horizontal line inside the box. The 90% and 10% percentiles are 
represented with vertical bars extending outside the box, and any outlying data points outside the 90% 
and 10% percentiles are individually represented by small circles. Thus, the relative ability of the eight 
failure models to predict failure was examined by comparing the predicted failure value in the different 
experimental tests for each of the materials, and by assessing if the predicted failure values for each 



material tended to have a constant value at the observed failure instances. 




To evaluate the different failure models, it is necessary to know the stress and strain 
distributions in the test specimens as a function of the applied deformation. For the uniaxial tension 
tests, the stress and strain distributions are directly obtained from the experimentally determined force- 
displacement curves. For the small punch tests, we obtained the stress-strain response by simulating 
the experimental tests using the finite element (FE) method. The FE simulations were based on the 
Hybrid Model (HM) [7-9] to represent the material behavior. 

3. Results 

The capabilities of the HM to predict the behavior of the UHMWPE is shown in Figures 2(a) 
and (b). These figures present a direct comparison between experimental data for the highly 
crosslinked UHMWPE material that was annealed at 1 10°C and predictions from the HM. Figure 3(a) 
shows the response in uniaxial tension at three different deformation rates. Figure 3(b) shows the 
force-displacement response from five repeated small punch test performed at a punch displacement 
rate of 0.5 mm/min. The failure points, indicated with the symbol □ in the figures, occurred at 
different applied deformations, but the experimental force-displacement responses are all very 
repeatable. As discussed by Bergstrom et al. [9-11] and illustrated in Figure 3, the HM is an accurate 
constitutive model for predicting the behavior of the different types of UHMWPE investigated in this 
work.. It can be seen that the failure strengths of the two highly crosslinked materials in uniaxial 
tension are lower than for the two conventional materials, but that the highly crosslinked materials fail 
at higher forces than the conventional materials when tested in small punch (equibiaxial) loading. This 
behavior is caused by different evolutions of the microstructural deformations when tested in uniaxial 
and multiaxial tests [11]. 

The results from applying the eight failure models to the experimental data for the four 
UHMWPE materials are summarized in Table 1. The results from applying the different failure 
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models to the experimental data are also graphically presented in box and whisker plots in Figures 4(a) 
to 4(h). The results for the maximum principal stress, Mises stress, and Tresca stress models show that 
these three failure models behave similarly (Table 1) and have almost identical failure conditions 
(Figures 4(a), (b), (c), respectively). They also show good predictive agreement between the uniaxial 
and the multiaxial failure data for all materials except the highly crosslinked material that was heat 
treated at 150°C. For this latter material, the effective stress value is significantly higher in the small 
punch tests than in the uniaxial tests. There is also a large scatter in the effective stress values at 
failure for these three failure models (Figures 4(a), (b), (c)). The variation in stress at failure for the 
unirradiated UHMWPE tested in uniaxial tension to failure, for example, varies between 150 and 320 
MPa. 

The predicted hydrostatic stress in the specimen at the time of failure is significantly higher in 
the small punch tests than in the uniaxial tests (Figure 4(d)). Hence, the hydrostatic stress model is not 
useful as a general failure predictor of UHMWPE. The Coulomb model is both qualitatively and 
quantitatively similar to the Mises, Tresca, and maximum principal stress models (Figure 4(e)). For all 
materials, the maximum principal strain at failure is larger in uniaxial tension than in the small punch 
experiments (Figure 4(f)). This is reasonable since the strain state at the tip of the small punch 
specimen at the time of failure is close to biaxial, and typically lower than in uniaxial tension. The 
experimentally determined Mises strain at failure is substantially different in uniaxial tension and small 
punch loading, clearly indicating that the Mises strain is not an accurate general failure predictor 
(Figure 4(g)). 

In contrast to the other failure models, using the chain stretch failure model, there is very good 
agreement between the predicted failure in uniaxial tension and in small punch loading, for all four 
materials (Figure 4(h)). Further, the variability in the critical chain stretch between different tests 
performed at different loading rates is small, relative to the other failure models (Figure 4(h)). The 
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value of the limiting chain stretch (2“), a material parameter in the HM that specifies how much the 
molecules can stretch before becoming fully extended is also shown in Figure 4(h). It appears that 
failure occurs when the molecular chain stretch reaches a certain fraction (60%-80%) of the fully 
extended chain stretch (Table 2). 

A direct comparison between the quantitative performances of the different failure models, by 
examination of the revalues between all experimental failure data and the critical failure values for the 
four materials demonstrates that the predictions from the chain stretch model have a much lower 
variability between different tests and better overall agreement with both uniaxial and multiaxial 
experimental failure data than the other failure models (Figure 5). The revalue for the chain stretch 
model is about 0.85. That is, the chain stretch model describes 85% of the variability in the failure 
data for a broad range of UHMWPE materials, loading conditions, and loading rates. Of the other 
models, the Mises stress, the maximum principal stress, Tresca stress, Coulomb stress and maximum 
principal strain are about equal, each having a r 2 -value of about 0.55. The hydrostatic stress and the 
Mises strain models are clearly the poorest performers, both having revalues less than 0.2. 

4. Discussion 

The experimental results and model predictions presented in this report indicate that uniaxial 
and multiaxial failure of conventional and highly crosslinked UHMWPE under monotonic loading can 
be accurately predicted using the chain stretch model. The chain stretch model is interesting not only 
because it provides the most accurate predictions of the models tested, but also due to its relationship 
to the stretching of the macromolecules in the material and with the maximum chain locking stretch 
parameter in the Hybrid model. The correlation between the experimental chain stretch at failure and 
the chain locking parameter suggests that if the HM is calibrated for a given UHMWPE material, then 
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a good multiaxial failure condition can likely be obtained by taking the chain stretch at failure to be 
about 60% to 80% of the chain locking stretch parameter. 

In practical terms, the chain stretch at failure for UHMWPE materials can be directly measured 
during standard uniaxial tensile tests. The chain stretch at failure is related to the uniaxial strain at 
failure (e) by the following equation: 




Uniaxial tensile tests are a standard method for characterizing the mechanical properties for 
conventional and highly crosslinked UHMWPE (e.g., ASTM F-648). Consequently, the chain stretch 
failure criterion outlined in the present study can readily be calculated for contemporary UHMWPE 
materials based on available experimental data. 

Some of the other failure models, e.g. the maximum principal stress and the Mises stress 
models, provide reasonable failure predictions of the uniaxial and small punch data. The scatter in the 
experimental failure criteria from these models, however, makes them less useful as a failure predictor 
under multiaxial loading conditions. 

The results from this study also suggest that the failure conditions for conventional and highly 
crosslinked UHMWPE depend on radiation dose and heat treatment (Table 1). One particularly 
interesting trend is observed when comparing the experimental failure data with the failure predictions 
of the different models. The failure of the unirradiated and the sterilized materials are both nicely 
captured by either a stress-based model (maximum principal stress, Mises stress, or Tresca stress) or 
by the chain stretch model. However, for the two highly crosslinked materials (particularly the 
remelted material 100 kGy, 150°C), the stress-based failure models tend to underestimate the strength 



12 




in the small punch experiments and overestimate the strength in the uniaxial experiments. The only 
model capable of accurately predicting the behavior of these materials is the chain stretch model. 

It is also interesting to consider the molecular weight between crosslinks and how it may relate 
to the maximum chain stretch. The concept of a molecular weight between crosslinks is derived from 
Flory’s [14] statistical theories for a swollen crosslinked polymer network in equilibrium with a good 
solvent. In contrast, the maximum chain stretch in our failure theory reflects the amount of local 
deformation that can be accommodated by UHMWPE, regardless of whether or not it is crosslinked. 
For virgin UHMWPE, the molecular weight between crosslinks is theoretically infinite, thus, there is 
no relationship between molecular weight between crosslinks and maximum chain stretch for the 
unirradiated UHMWPE examined in this study. With respect to the materials considered in our study, 
the molecular weight between crosslinks decreases from infinity (for unirradiated UHMWPE) to an 
estimated value of 6,000 g/mol at a dose of 100 kGy [28]. For this reason, we do not judge the 
molecular weight between crosslinks to be an attractive state variable for predicting the failure 
behavior of conventional and highly crosslinked UHMWPE. 

It is known that the material response of UHMWPE is rate-dependent [9]. By performing an 
analysis of variance (ANOVA) on the experimental data obtained at different loading rates it is 
possible to examine if the failure quantities are rate-dependent. A direct analysis of the experimental 
data shows that the only rate-effects that are statistically significant (p < 0.05) are: (1) the punch force 
at failure increases with increasing punch rate; (2) the punch displacement at failure decreases with 
increasing punch rate for the unirradiated material. The other failure quantities, e.g. the stress and 
strain response in uniaxial tension, showed no statistically significant rate-dependence [20]. Overall, 
these findings support our assumption that the failure criteria for UHMWPE is effectively rate 
independent. 
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The high stress and strain values at failure that are reported in Table 1 are indicative of the 
intrinsic strength of the material and can be used to predict material failure on the microscale at the tip 
of microcracks or the local damage during wear. A complicating factor when studying failure of 
polymeric components is that failure of seemingly identical components typically occurs at different 
load levels. The apparent distribution in failure strengths is caused by variations in the internal state of 
the material in different components. The differences can, for example, be due to localized differences 
in molecular weight, molecular weight distribution, crosslinking density, degree of crystallinity, 
morphology, or defect population. The failure on a macroscopic scale therefore often appears 
stochastic. The results presented here for monotonic loading lay the foundation for the future 
development of more advanced stochastic damage models that can accurately predict failure under 
multiaxial, cyclic loading conditions. 
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Appendix A. 

This appendix briefly summarizes the governing equations of the Hybrid Model (HM). For a 
more in depth discussion, see the recent publications by Bergstrom et al. [9, 10]. 

The kinematic framework used in the HM is based on a conceptual homogenization of the 
semi-crystalline and amorphous domains [2, 3, 16], and a decomposition of the deformation gradient 
into elastic and viscoplastic components: F = F'’ F p . The Cauchy stress acting on the system can be 
obtained from the linear elastic relationship: 

T + A, tr[E c ]l), (al) 
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where E e = ln[V'J is the logarithmic true strain, f = det[F e ], and [f, X e are Lame’s constants. For 
convenience, the Lame constants can be replaced by the equivalent Young’s modulus (E) and 
Poisson’s ratio ( v):ju e = E e l (2(1 +v e )) and X e = E e v e / ((l+v e )(l-2v e )). 

The details of the stress expressions are presented in Equations (a2) - (a4) 

T, = T 2^[T w „ i ,,(F'';^,Ar,*'d+<?A(F'-; lUj )], (a2) 



T = 

a 8 chain 



r 



L l 

^(i >K ck ) 



dev]^], 



(a3) 



(a4) 

where Ta is a tensor-valued function of the viscoplastic deformation gradient F p and the material 
parameters (pa, XA lock , ka, ^a}, where pa is the shear modulus, XA lock is the locking stretch, k a is the bulk 
modulus, and qA is a material parameter specifying the relative magnitudes of T schain and T I2 , and B p * is 
the left Cauchy-Green deformation tensor, 



The stress driving the viscoplastic flow of the backstress network is obtained from the same 
hyperelastic representation that was used to calculate the backstress: T B = s B • T, ) , 

where s B is a dimensionless material parameter specifying the relative stiffness of the backstress 
network. The parameter s B evolves with imposed plastic deformation to capture the distributed 
yielding: s B = p B - (,s fi - s Bf where p B is a material parameter specifying the transition rate of the 

distributed yielding event, s B f is the final value of s B reached at fully developed plastic flow, and f' B is 
the rate of viscoplastic flow of the time-dependent part of network B: 
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f B =to 



(a5) 



v " b J 

where Z B = ||dev [T b ]||^ , z b “ se and m B are material parameters, and j 0 is a constant coefficient with a 
value of 1/s. The velocity gradient of the viscoelastic flow of the backstress network is given by 






(a6) 



The yielding and plastic flow of the material is captured using: 



L P =fc 



dev[T c ] R / 
T c ] 



(a7) 



where If - F P F P ' , T c = T -[F e (T A + T B )F eT ]/ J e is the stress acting on the relaxed configuration 
convected to the current configuration, T c — |dev[T r j| f . is the effective shear stress (calculated using 
the Frobenius norm) driving the viscoplastic flow, f c = f Q • {z c / T^ ase ) m<:: is the magnitude of the 
viscoplastic flow, z basr and me are material parameters, and f 0 a constant coefficient with a value of 
1/s. 

In total, the augmented HM contains 13 material parameters: 2 small strain elastic constants ( E e , v e ); 4 
hyperelastic constants for the back stress network (p A , ^A° ck , k a , qf)\ 5 flow constants of the backstress 
network (s B i, s B f, p B , z B base , m B ); and 2 yield and viscoplastic flow parameters (xc base , me). Of these 13 
material parameters, four (elastic modulus (£); yield strength ( Z b B ase ); the effective stiffness after yield 
( p A ); and the limiting chain stretch ( X '"‘ :k ), which controls the large strain behavior) have been shown 
to vary between UHMWPE formulations, while the remaining nine are constant. The material 
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constants used in the analysis are listed in Table 3. The HM as presented here has, to date, been 
validated for UHMWPE formulations within the crystallinity range of 46% to 61% and for both 
uncrosslinked and crosslinked (100 kGy Y-N 2 ) UHMWPE. 
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Figure 2. (a) Uniaxial tension to failure data, (b) small punch loading to failure. The experimental tests were 
performed on UHMWPE prepared by four different thermal and radiation treatments. 
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Figure 3. Comparison between the experimental data and Hybrid Model predictions: (a) uniaxial tension to 
failure, (b) small punch testing to failure. The experiments were performed on highly crosslinked 
UHMWPE annealed at 1 10°C. The failure points are indicated by the symbol □. 
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Figure 4. Predicted failure conditions in uniaxial tension and small punch loading of the four different types of 



UHMWPE examined: (a) maximum principal stress; (b) Mises stress; (c) Tresca stress; (d) hydrostatic 
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Table 1. Experimentally determined failure conditions for the eight failure models and the four types of 



UHMWPE. 



Failure 

Condition 


Unirradiated 


30 kGy 
(Y-N 2 ) 


100 kGy 

(110°C) 


100 kGy 

(150°C) 


Max principal 


246.2 


204.1 


172.1 


158.2 


stress (MPa) 


(39.7) 


(25.6) 


(25.4) 


(33.8) 


Mises stress 


242.1 


200.6 


170.3 


156.6 


(MPa) 


(39.1) 


(25.6) 


(24.4) 


(31.3) 


Tresca stress 


246.9 


204.4 


172.6 


158.6 


(MPa) 


(39.9) 


(25.7) 


(25.6) 


(34.0) 


Hydrostatic stress 


120.8 


99.2 


85.9 


81.1 


(MPa) 


(44.8) 


(33.5) 


(35.4) 


(39.6) 


Coulomb stress 


271.8 


224.8 


190.1 


174.4 


(MPa) 


(44.3) 


(28.2) 


(28.6) 


(37.4) 


Max principal 


1.48 


1.37 


1.12 


1.13 


strain 


(0.18) 


(0.19) 


(0.14) 


(0.09) 


Mises strain 


1.61 


1.46 


1.12 


1.19 


(0.80) 


(0.70) 


(0.60) 


(0.59) 


Chain stretch 


2.90 

(0.22) 


2.60 

(0.21) 


2.04 

(0.08) 


2.01 

(0.08) 
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Table 2. Comparison between the experimentally observed chain stretch at failure and the chain locking stretch 



parameter. 





Unirradiated 


30 kGy 
(Y-N 2 ) 


100 kGy 
(110°C) 


100 kGy 
(150°C) 


Chain stretch 
at failure 


2.90 


2.60 


2.04 


2.01 


Chain locking 
stretch 


4.44 


4.40 


2.80 


2.52 


parameter 
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Table 3. Hybrid model material parameters for each of the four UHMWPE materials used in the 



tension and small punch simulations. 



Material 

Parameter 


Unirradiated 


30 kGy 

y-n 2 


100 kGy 
110°C 


100 kGy 
150°C 


E e (MPa) 


1574 


2020 


2009 


1270 


V e 


0.46 


0.46 


0.46 


0.46 


ptA (MPa) 


7.73 


8.22 


10.15 


8.14 


Mf 


4.44 


4.40 


2.80 


2.52 


k a (MPa) 


2000 


2000 


2000 


2000 


qA 


0.20 


0.20 


0.20 


0.20 


SBi 


40.0 


40.0 


40.0 


40.0 


S B f 


10.0 


10.0 


10.0 


10.0 


Pb 


27.0 


27.0 


27.0 


27.0 


T b B ase (MPa) 


26.1 


25.0 


26.2 


20.7 


m B 


9.5 


9.5 


9.5 


9.5 


T b c ase (MPa) 


8.0 


8.0 


8.5 


8.5 


me 


3.3 


3.3 


3.3 


3.3 
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